A SIMPLE COMPUTATION OF C(2fc) BY USING BERNOULLI 
POLYNOMIALS AND A TELESCOPING SERIES 
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Abstract. We present a new proof of Eulcr's formulas for C(2fc), where k = 
1, 2, 3, ... , which uses only the defining properties of the Bernoulli polyno- 
mials, obtaining the value of C(2fc) by summing a telescoping series. Only 
basic techniques from Calculus are needed to carry out the computation. The 
method also applies to C{2k + 1) and the harmonic numbers, yielding integral 
formulas for these. 



1. Introduction 

In the mathematical literature, one finds many ways of obtaining the formula 



(1) Cm:=J2^= ' " (2fc)! " fc = 1,2,3,..., 

n—l ^ ' 

where Sfc is the fcth Bernoulli number, a result first published by Euler in 1734. 
For example, the recent paper [1] contains quite a complete list of references. The 
aim of this paper is to give a new proof of (1) which is simple and elementary, in 
the sense that it involves only the Bernoulli polynomials and a telescoping series. 

2. The Bernoulli polynomials 

For completeness, we begin by recalling the definition of the Bernoulli polyno- 
mials B}^(t) and their basic properties. There are of course multiple approaches 
one can take (see [4], which shows seven ways of defining these polynomials). A 
frequent starting point is the generating function 



1 ^_2^^fc-l22A;-1^2fe 



(2) ;^-E^^w 



k=0 

from which, by the uniqueness theorem for power series expansions, one can quickly 
obtain many of the basic properties, among which we single out 

(3) 5oW = l, B'kit) ^ kBk^iit), fc>l, 

which by induction shows that Bk{t) is in fact a polynomial, and 



(4) I Bk{t)dt = 0, k>l. 

Jo 

If one considers that the generating function approach is not completely "elemen- 
tary," in the sense of requiring more than standard first year Calculus (series of 
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functions, justifying exchanging the derivative with the infinite sum, etc.) it is eas- 
ily avoided by instead using (3) and (4) to define the polynomials Bk{t). In any 
case, one finds that the first few Bernoulli polynomials are 

The Bernoulli numbers are defined to be the values Bk = Bk{0). From (3), (4), 
and the Fundamental Theorem of Calculus, one sees that Bk{Q) = Bk{l) for k > 2, 
and from the symmetry relation 

(5) -t) = (-l)^Sfe(t), fc>0 

(easily proved by induction on k) one deduces 

B2kil) = B2ki0) = B2k, k>0, 

B2k+i{l) = B2k+iiO)=0, k>l. 

There are many other properties and relations satisfied by the Bernoulli polyno- 
mials, but those listed above are the only ones necessary for the goal we have in 
mind. 



3. Computation of C{2k) 
3.1. Some auxiliary integrals. Consider the integrals 

I{k,m) := / B2kit) cos{rm:t) dt, k > 0, m > 1. 
Jo 

An immediate computation shows that 1(0, m) ~ for to > 1, and, for k > 1, 
integrating by parts twice and applying (3), we get 

1 r 2k 
I(k,m,) = B2k(t) sm(mTTt) / B2k-iit) sm(imTt) dt 

TOTT L J t=0 TOTT Jo 



2k 



TO^TT^ 

which gives us both the special case 



=1 2k(2k — 1) 

B2k-lit) COs{mTTt) — /(fc-l,TO), 




TO^TT 



/•I/ i\ |0, m = l,3,5,. 

1(1, to) = / [ t" — t -\ — cos(TO7rt) dt 
Jo \ GJ 

and, by (6), the recurrence relation 

/(fc,TO) = -^^^^J^/(fc-l,TO), k>2, 

from which, by induction, one obtains the closed form 

{0, m=l,3,5,..., 
91, 91, 7 J/t — ^, 1, u, . . . . 
TO TT 

Now, for reasons which are made clear below, consider B'^{t) = Bk{t) — Bk{Q) = 
Bk{t) — Bk, i.e., the Bernoulli polynomial minus its constant term. The correspond- 
ing integral 

r{k, to) ;= / B2k{t) cos(TO7rt) dt = [ [B2k{t) - B2k) cos{rmrt) dt 
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is equal to I{k, m), because cos{m'Kt) dt = 0. For fixed k > I, summing I*{k, m) 
over m > 1 yields, by (7), 



rn— 1 



Tn=l ^ ' 

oo oo oo 

= ^ r(fc,2m) + ^ r(fc,2?72 + 1) = ^ I*{k,m) 



rn — 1 Tn — ?n — 1 

oo ^1 



/ B2fcWcos(m^t)dt. 

rn=l 



3.2. The telescoping trick. We need the elementary trigonometric identity 

^ ^ sini^^x) - sin(2H|:^x) 
(8) cos(m.) = ^—^^ , 

which can be proved, for example, by induction. With the introduction of (8), we 
now have a telescoping scries, yielding 

^^°°r^i 2sm( — j Jo ^sm( — ) J 

We observe that by (4), the value of the second term is 

\ J^' B;,it) dt^\ W - S2.) dt = 

Now, we show that the limit in the first term is zero. Note that the function 

extends by continuity to t = since i32/c(0) (this is the reason for subtracting 
the constant term), and is differentiable on [0,1] with a continuous derivative. 
Denoting (2A^ + l)7r/2 by R, integrating by parts gives 

m sHRt) dt . -2^^/(1) + 1/(0) + 1 fitf-^ dt. 

The boundedness of f'{t) and the trivial estimate | cost| < 1 show that each term 
in the above sum approaches zero as i? ^ oo, so that indeed the limit tends to 0. 
Consequently, 

22fc^2fc U^fc)- 2 
which after rearranging terms, gives (1). 
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4. Fourier Confidential 

Readers familiar with the basic theory of Fourier series will of course immediately 
recognize that the integrals /(fc,m) are essentially the Fourier coefficients of the 
Bernoulli polynomial i?2fe(0- Namely, the Fourier series of B2k {t), first computed 
by Hurwitz, is 

and setting t = yields (1). Of course, this approach necessitates proving the 
basic facts about Fourier series, not the least of which are the issue of pointwise 
convergence and the Riemann-Lebesgue Lemma. Fortunately, the latter is easily 
proved for functions with the argument using integration by parts we have given, 
and the issue of convergence is dealt with using the "telescoping trick" (8). 

The use of the same "telescope" (8) in [2] brought to our attention that this 
idea, used there to compute C(2) via the integral xcos{mx)dx, could also be 
used to compute C{2k). However, one needs to realize that instead of integrating 
x'^cos(mx), the natural functions to integrate against the cosine are the Bernoulli 
polynomials. 

The difference formula (8) will be recognized as that obtained from the definition 
of the Dirichlet kernel and its expression as 

^ ^ sin f ^AT+i ^j-) 

DN{t)= Yl * = 1 + 2 ^ cos(m7rO = ^ ' > . 

1 sm ^ 

m=—N m=l \ Z / 

It is also worthwhile noting that a related telescoping trick also yields a nice result 
regarding pointwise convergence of Fourier series in general, as shown in [3]. 

One also needs to see how the basic defining relations such as (3) can be employed 
to give an elementary method of computation for the Fourier coefficients. We have 
used the recurrence relation (3) to derive the formulas for the coefficients in a 
straightforward manner. This approach is actually not standard, although it seems 
to us to be the simplest. Most classical approaches involve the use of complex 
analytic summation techniques, due in part to the fact that they generalize, for 
example, to the Hurwitz zeta function. Such techniques are quite useful in analytic 
number theory, but would negate the elementary nature of the proof given here. 

5. What about Q{2k + 1)? 

The same approach will yield a formula for Q{2k + 1), but the term which is 
subtracted when summing the telescoping series is an integral which, as far as 
anyone knows, cannot be evaluated in (a non-tautological) closed form. We proceed 
in the same way as before, except this time, we consider the integrals 

J(fc,m):= / i32fe+i(t) sin(TO7r<) o?i 
(which are the Fourier sine coefficients). Direct computation shows that 

/•I/ i\ 1 + (-!)'" fo. m = 1,3,5,..., 

J(0, m) = / t - - sinfmTTi) dt = i '— = { 1 
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and, for k > 1, integrating by parts twice gives, using (3) and (6), the recurrence 
relation 

(2fc + l)(2fc) 
Jffc, m) = J{k - 1, mj 

from which, by induction, we obtain the closed form 

0, m = 1, 3, 5, . . . , 

(9) J(/c,m) ^ ( (-l)fc-i(2fc + l)! 

' -J — , m ~ 2,4,6, 

Note that this time, since i?2fe+i(0) ~ -B2fe+i = for A: > 1, subtracting the constant 
term is not necessary except for fc = 0, since Bi = —1/2, although, since C(l) = oo, 
this is irrelevant for now (but see the next section). Thus, summing (9) over m, we 
obtain, for fc > 1, 



oo 



m— 1 m—l 

oo oo oo 

= ^ J{k, 2m) + ^ J{k, 2m + 1) = ^ J{k, m) 

m— 1 ?n— m—l 

oo „i 

= / i32fc+i(t) sin(m7rt) di. 

The trigonometric identity that plays the role of (8) in this case is 

, , cos(2h^x) - cos(22^a;) 
(10) sin(m.) = - - ^ ^ 

Thus, we now obtain the telescoping series 

/ C0s(2im^t) \ /•! , C0s(f ) 

lim / B2k+i{t) , . ' rft + / g2fc+i(t) , . 7,;, dt. 



\^iv^oo7o --r-v / 2sin(f) j Jo ^ "2sm{f) 
The limit is null, for the same reason as before. Consequently, we have proved that 

f-l^'=-i22'=7r2'=+i /Trt\ 
(11) C(2fc + 1)= ^ {2k + l)\ i ^2fe+i(t)cot(-)dt, k>l. 

For example, 

.3 /•! 



C(3) = ^^ (2t3-3i2+t)cot(|)dt. 



It would be nice to know if the integral in (11) has a closed form expression (other 
than in terms of C(2fc + 1) of course!) but at present this problem remains open. 
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6. C(l): THE HARMONIC SERIES 

Although the harmonic series diverges, we can still obtain information from the 
integrals J(0, m) if we consider the partial sums of the telescoping series in the last 
section with fc = 0. It is actually easier not to remove the constant term Bi = —1/2, 
but if we don't, note that each difference of integrals in the telescoping sum is now 
an indeterminacy of the form oo — cxd. The trick is to just use (10) to sum sin(m7ra;) 
inside the integral J(0,to). Then it is easy to see that we obtain the formula 

/"V 1\ cos(Mm^i) _cos(f) 1 
Jo V 2y sm(^) TT 

where Hm = Y^^lf , — is the Mth harmonic number. 

By (11), the same formula is valid with AM + 3 instead of AM + 1, but for AM 
and AM + 2, the interested reader may verify that we get combinations involving tt, 
log 2, and Catalan's constant 

(-1)" 



G 



OO 

^ (2?i+ 1)2 



via the formulas 



^ t , 8 ^ /ttA , 21og2 
■dt^ ^G, / t cotl — ] dt^ 



/o sin(f) 7r2 ' J, \2 

but these, though certainly interesting, are not really part of this story. 
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